UNCLASSIFIED 


AD  NUMBER 


AD482327 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  Feb  1966. 
Other  requests  shall  be  referred  to  Army 
Electronics  Command,  Atmospheric  Sciences 
Lab,  White  Sands  Missile  Range,  NM 
88002-5501. 

AUTHORITY 


US  Army  Electronics  Command  ltr  dtd  20  Sep 
1971 


THIS  PAGE  IS  UNCLASSIFIED 


DDC  AVAILABILITY 


Qualified  requesters  may  obtain  copies  of  this  report  from  DDC. 


DISPOSITION  INSTRUCTIONS 


Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it  to  the  originator. 


DISCLAIMER 


The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of  the 
ny  position,  unless  bo  designated  by  other  authorized  documents. 


CONTINUOUS  DETERMINATION  CF  THE  AVERAGE  SOUND 


(W) 

DA  Tft6K^/jP620901Al98-02 


ATMOSPHERIC  SCIENCES  LABORATORY  <j0Q 
WHITE  SANDS  MISSILE  RANGE,  NEW  MEXICO 


ABSTRACT 


A  reasonably  accurate  determination  ot  the  velocity  of  sound  In  the  open 
air  along  a  path  of  arbitrary  length  by  some  inexpensive  means  Is  the  crux  of 
sonic  atmospheric  sensing.  At  present,  this  is  accomplished  by  measuring  either 
time  delays  or  phase  differences  of  an  acoustic  signal.  This  paper  explores  the 
technique  of  phase  measurement  and  proposes  a  direct  approach  to  the  problem  of 
determining  the  sound  speed  within  an  arbitrary  layer.  By  means  of  the  phase 
technique,  the  initial  conditions  of  the  sound  velocity  and  the  vector  wind  may 
be  determined.  The  temporal  distribution  of  the  sound  speed  is  found  from  an 
equation  for  the  average  sound  velocity  along  the  paths  between  the  sound  source 
and  two  opposing  transducers. 
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INTRODUCTION 


The  fundamental  problem  in  tonic  atmospheric  tenting  lies  in  the  temporal 
determination  of  an  accurate  speed  of  sound  in  the  open  air  and  the  separation 
of  the  wind  effects  from  the  scalar  thermodynamic  sound  speed.  It  is  apparent 
that  the  accuracy  of  the  latter  is  highly  dependent  upon  that  of  the  former. 

The  usual  technique  (Barrett  and  Suomi,  1949;  and  John,  i960)  is  to  measure  the 
difference  and  sum  of  transit  times  of  the  sound  pulses  from  a  source  travelling 
in  opposite  directions  to  two  transducers,  each  at  a  specified  distance  from 
the  source.  At  high  frequencies  and  short  distances,  the  difference  in  transit 
times  Is  generally  of  the  order  of  microseconds.  However,  the  same  result  may 
be  accomplished  by  measuring  the  relative  phase  angles  (Coffman  and  Price,  1962), 
thereby  magnifying  the  sum  and  difference  by  a  factor  of  2  f,  .heir  f  is  the 
frequency  used.  The  phase  technique  was  utilized  over  a  long  path,  approximate¬ 
ly  60  meters  (Low,  1964),  and  was  found  to  give  satisfactory  results  when  com¬ 
pared  to  the  results  derived  from  an  array  of  thermocouples  along  the  same  path. 
Kaimal  and  Businger's  system  (1963),  constructed  on  the  same  principle,  measures 
the  vertical  wind  and  temperature  variations  approximately  at  a  point  by  artifi¬ 
cially  presetting  certain  mean  values  of  the  temperature  and  hence  of  the  sound 
velocity  under  windless  conditions. 

In  the  present  paper,  the  technique  of  phase  measurement  in  sonic  atmos¬ 
pheric  sensing  is  extended.  Instead  of  attempting  to  obtain  solutions  for  the 
wind  components  from  the  basic  equations,  an  effort  is  made  to  determine  con¬ 
tinuously  a  reasonably  accurate  value  of  the  sound  velocity  in  the  open  air  over 
an  arbitrary  path.  Then  it  is  a  simple  matter  to  extract  both  the  wind  and 
temperature  data  along  the  same  path. 

THEORY 

Consider  the  isoscles  triangle  ABC  in  Figure  1  and  assume  that  the  wind 
vector  can  be  resolved  along  the  path  of  sound  propagation,  bearing  in  mind 
Schotland's  (1955)  excellent  analysis  of  wind  measurement  by  sonic  means.  A 
sound  source  is  placed  at  the  vertex  A  and  a  transducer  at  each  of  points  B  and 
C.  In  addition  to  those  labeled  in  the  figure,  other  notations  are  listed  below; 


Figure  1 
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D  =  AB  =  AC 

C0,  C  =  the  initial  and  subsequent  velocities  of  sound; 

{,  {'  -  the  frequencies  used  at  the  sound  souiee; 

in,  n  =  the  numbers  of  wavelengths  in  the  paths  AB  and  AC  respectively 

at  the  given  frequencies; 

t  ,  t  =  the  initial  and  subsequent  transit  limes  from  A  «.o  B; 

t^,  t'  =  the  initial  and  subsequent  transit  times  from  A  to  C; 

M0l  N  =  the  initial  and  subsequent  wind  components; 

>.,  A'  sc  the  initial  wavelengths  along  the  paths  AB  and  AC  respectively; 

£♦,  =  the  phase  differences  along  the  paths  AB  and  AC  respectively. 

Symmetry  has  been  assumed  (see  Figure  1  and  notations). 

The  initial  conditions  require  that  for  the  frequencies  f  and  f*  used,  the 
numbers  m  and  n  must  be  accurately  determined  such  that: 


for 

path 

03 

D  =  lex 

and 

(1) 

f  A  = 

C0  -W0  cosa  } 

(2) 

for 

path 

AC 

D  s  n*' 

and 

(3) 

f  1 A  1 

=  C  +  W0  cosa. 

(4) 

Under  these  conditions,  it  is  preferable,  though  not  necessary,  to  make  m 
and  n  integers. 

The  travel  time  of  a  wave  front  emitted  from  the  sound  source  A  to  point  B, 
by  utilizing  equations  (1)  and  (2),  is  given  by 

tQ  =  D/(C0  -  V0  coaa).  (5) 

The  sound  velocity  and  the  wind  component  are  now  denoted  by  C  and  W,  re¬ 
spectively,  and  since  they  no  longer  represent  initial  conditions,  a  new  travel 
time  is  obtained 

t  =  D/(C  -  W  co ao).  (6) 


The  difference,  At  =  t0  -  t,  is  the  time  delay  or  time  difference.  From 
equations  (5)  and  (6)  the  time  difference  may  be  determined  as 

At=  - £ -  -  - E - -  (7) 

C0  -  Wc  coaa  C-  V  eoao 

Solving  for  the  velocity  C  by  a  series  of  algebraic  manipulations,  it  is 
found  that 
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♦  W  COM2. 


C  -  W  com 

o  o 

(c„  -w  com  la  t 

,  o  O 


Equation  (6)  nay  be  simplified  by  means  of  equations  (1)  and  (2),  noting 
that  f/m  =  (C0  -  h'0  cos  a  )/D, 


ooea 


'+  W  coma- 


1  -  £At/« 


<9) 


The  difference  in  time  At  can  now  be  expressed  as  the  difference  in  phase 
angle  A+  by  the  formula  A+  =  2  n  f  At.  Substitution  of  this  formula  into 
equation  (9)  gives  the  desired  expression  for  the  average  sound  velocity  C  in 
terms  of  phase  difference  with  respect  to  the  initial  conditions 


C 


C  -W  coma 

O  O 

1-  A#7/aa  n 


♦  V  coaa* 


(10) 


By  the  same  series  of  manipulations,  using  At'  =  t^  -  t'  for  the  path  AC. 
the  following  expression  for  the  average  sound  velocity  along  AC  is  found  to  be: 


C  +  W  coaa 
C  =  o  ° 

1  -  a*’  /2nn 


U  cos  a  . 


(ID 


Expressions  (10)  and  (11)  form  a  oair  of  simultaneous  equations  and  can  be 
solved  exactly,  if  so  desired;  however,  the  result  may  not  sufficiently  justify 
the  effort.  Solving  for  the  average  sound  velocity  C  by  combining  the  above  two 
equations,  we  obtain 


2C  -  (Co  -Vo  ooaa)/(l  - 

♦  (Co  ♦  Wo  eo*a)/(l  -  A*'/2nI  ).  (12) 

I 

Generally,  A* ,  A#  <«2ar  ,  2n«'  respectively;  it  is  reasonable 

to  take  the  following  approximations: 

1/(1  -  -  1  A»/2n*,  and 

1/(1  -  A*'/2nl)  •  1  ♦  A»'/2nI*  (]3) 

Upon  insion  of  equation  (12)  using  the  above  approximations  and  factor¬ 
ing,  we  o.  -in  two  terms 
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-(A*/2ml)  W  coso 4\nd  4-  (M'/Snl)  W  coaQ^C 

'  O  O  O  O 

(14) 

which  are  not  much  different  from  zero.  Thus,  a  reasonably  accurate  expression 
for  the  continuous  average  sound  velocity  along  an  arbitrary  path  is  given  by 

C  -  Co  [1  +  (A»/m  +  A*7n)].  (15) 

The  average  acoustic  virtual  temperature  (Suomi,  1949)  in  degrees  absolute 
along  the  path,  which  will  not  be  appreciably  different  from  the  average  free 
air  temperature  along  the  same  path  with  low  humidity  conditions,  is  readily 
obtainable  from  the  average  sound  velocity  C  given  above  by  means  of  the  Laplace 
formula . 

To  obtain  the  average  horizontal  wind  component  across  the  paths,  equations 
(10)  and  (11)  are  solved  simultaneously  and,  with  approximations,  yield 

ZW  cosa  «=  2  V  cosa  -  (C  -  W  eoea) 

C  0  0 

+  (C  +  W  cosak*  '/2nB> 
o  o 


(16) 

Equati  ons  (1),  (2),  (3),  and  (4),  the  initial  conditions,  and  the  fact  that 
cosa  =  d/"D  are  used  to  simplify  this  expression.  The  average  horizontal  wind 
component  h'  is  then  given  by 

w  ;  wo  -  isr  n7) 

DISCUSSION  AND  CONCLUSIONS 

The  initial  conditions  may  conceivably  cause  sume  concern  as  to  their  mean¬ 
ing.  Actually,  the  setting  of  the  initial  conditions  could  be  termed  calibra¬ 
tion  which  is  a  required  step  for  any  instrument.  Under  normal  weather  condi¬ 
tions,  this  setting  need  not  be  at  all  altered  if  there  is  no  drift  in  the  elec¬ 
tronic  equipment , 

In  the  course  of  developing  the  equations  in  the  preceding  section,  no  re¬ 
strictions  hav'  been  imposed  thereon,  except  for  approxinat  ions  which  are  con¬ 
sidered  reasonable  and  should  in  no  manner  jeoprrdize  the  precision  of  measure¬ 
ment  within  the  limitations  of  the  electronic  equipment.  Equations  (15)  and 
(17)  are  generalized  forms;  however,  the  frequencies  chosen  do  not  have  to  be 
different,  in  which  case  either  m  or  n  may  take  on  the  value  of  a  noninteger. 

The  use  of  a  single  frequency  does  not  greatly  simplify  the  cpmputations ,  but  it 
has  the  advantage  of  presenting  a  visual  picture  of  the  wind  effect  when  the 
temperature  remains  unchanged  during  a  short  interval  of  observation,  and  it  may 
bring  about  some  economy  in  instrumentation.  Furthermore,  as  nay  be  obvious 
from  equations  (2)  and  (4),  the  wind  effect  will  cause  the  phase  angles  and 
£♦’  to  assume  opposite  signs,  thereby  indicating  the  sense  of  the  wind  direc¬ 
tion,  while  the  temperature  effect  will  cause  them  to  have  the  same  sign,  either 
negative  or  positive,  depending  on  the  initial  calibration  values. 
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The  system  is  highly  sensitive  to  wind  and  temperature  changes;  '^erefore, 
the  frequency  or  frequencies  I"  b«  used  should  be  selected  with  <*- ..e,  Unless  a 
higher  resolution  is  required,  the  lower  audible  frequenc f».i  should  be  used  for 
reasons  noted  below.  As  shown  before,  &♦"  2  frit,  hrormil}  a  change  of  1°C  in 
temperature  will  cause  a  change  of  about  0.t>  meter  per  second  in  sound  speed. 

Now  assume  that  the  time  difference  At,  as  a  result  of  a  change  of  1°C,  is  nbout 
1  millisecond.  It  Is  easily  seen  that  it  f  -  1 00  cps,  =  0.36°;  at  f  *  lfiCO 
cps,  -  3.6°,  etc.  The  same  reasoning  holds  for  a  change  in  the  wind  field. 
Moreover,  the  longer  the  path  D,  the  larger  the  time  difference  and  the 
greater  the  desirability  to  use  lower  frequencies. 

As  the  angle  a  in  Figure  1  becomes  larger  and  larger,  the  value  of  W  coso 
approaches  the  vertical  component  of  the  vector  wind.  As  the  angle  o  approaches 
90  degrees,  the  system  will  have  become  a  device  for  temperature  measurement, 
the  vertical  components  being  negligible  under  normal  conditions.  The  average 
free  air  temperature  in  degrees  absolute  along  the  path  may  then  be  obtained 
directly  from  the  observed  phase  angle 


T  =  T0  (1  ♦  0.00556  A*/a). 


(10) 


At  the  other  extreme,  as  the  value  of  cosa  approaches  one,  the  system  as¬ 
sumes  the  familiar  geometry  employed  in  sonic  anemometry-thermometry .  As  the 
distance  between  the  transducers  and  sound  source  is  reauced,  approximate  point 
measurement  is  attained. 


It  may  be  recalled  that  in  the  text  the  word  "wind  component"  has  been  em¬ 
ployed  to  emphasize  the  fact  that  the  "Ind  so  obtained  is  not  the  total  wind 
vector  but  a  component  thereof,  depending  on  the  orienialiun  of  the  system.  If 
two  more  transducers  are  available  so  that  each  is  oriented  orthogonally  to  the 
next,  it  is  then  possible  to  obtain  the  average  total  wind  field  across  the 
array. 
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